■’^’Calhoun 

IflitiKuliarttl  Arehivtaf  (he  Naval  ftutgrdduH*  School 


Authors) 

Griggs,  J ohn C.;  Macdonald,  Alexander  R. 

Title 

Direct  decay  half-life  measurements  of  isometric  states  of  HF(178),  HF(179),  and 
YB(177) 

Publisher 

Monterey,  California:  U.S.  Naval  Postgraduate  School 

Issue  Date 

1963 

URL 

http://hdl .  handl  e  net/10945/12018 

This  document  was  downloaded  on  May  12,  2015  at  03:46:49 


DUDLEY 

KNOX 

LIBRARY 


htt  p:/ /www.  n  ps.  e  d  u/l  ib  ra  ry 


Calhoun  is  a  project  of  the  Dudley  Knox  library  at  NP5,  furthering  the  precepts  and 
goals  of  open  government:  and  government  transparency.  All  information  contained 
herein  has  been  approved  for  release  by  the  NPS  Public  Affairs  Officer. 

Dudley  Knox  Library  /  Naval  Postgraduate  School 
411  Dyer  Road  /  1  University  Circle 
Monterey,  California  USA  93943 


http://www.n  ps .  edu/ 


NPS  ARCHIVE 
1963 

GRIGGS,  3. 


DIRECT  DECAY  HAIF-UFE 
MEASUREMENTS  OF  ISOMERIC  STATES 
OF  HP™,  HFmt  AND  YB177 
JOHN  C.  SR!®©$ 

ALEXANDER  R.  MACDONALD 


DUDLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93943-5101 


$0 


DIRECT  DECAY  HALF-LIFE  MEASUREMENTS  OF 
ISOMERIC  STATES  OF  HF17^,  HF179S  AND  YB177 


#  #  #  # 

John  C.  Griggs 
and 


Alexander  R.  Macdonald 


DIRECT  DECAY  HALF-LIFE  MEASUREMENTS  OF 


ISOMERIC  STATES  OF  HF1^,  HF179j  AND  yb177 


by 

John  C.  Griggs 

/  * 

Captain,  United  States  Army 
and 

Alexander  R.  Macdonald 
Captain,  United  States  Army 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 

PHYSICS 

United  States  Naval  Postgraduate  School 
Monterey,  California 


19  6  3 


DUDLEY  KNOX  LIBRARY 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93943-5101 


DIRECT  DECAY  HALF-LIFE  MEASUREMENTS  OF 
ISOMERIC  STATES  OF  HF1?^  hfW5  AND  YB1?? 

by 

John  C„  Griggs 
and 

Alexander  R.  Macdonald 


This  work  is  accepted  as  fulfilling 
the  thesis  requirements  for  the  degree  of 
MASTER  OF  SCIENCE 
IN 

PHYSICS 
from  the 

United  States  Naval  Postgraduate  School 


ABSTRACT 


Half-lives  of  thermal  neutron  activated  samples  were 
determined  by  direct  decay  measurements,  Transport  of 
samples  between  activation  and  counting  areas  was  accomplished 
by  employment  of  a  rapid  pneumatic  transfer  system,  Use  of  a 
single  channel  analyser  with  a  fast  chart  recorder  to  record 
scaled  counts  vs.  time,  and  analysis  of  these  data  in  a 
computer  program  are  discussed.  Results  for  nuclides 
investigated : 

Hf17&n  4  s  (3<>92  t  .05)  sec0 

Hf179m  ti  .  (1&049  t  0O4)  sec, 

Yb177m  4  s  (6,49  t  0O3)  sec. 

Our  method  has  several  advantages  over  multichannel 
devices  and  appears  to  be  capable  of  measuring  half-lives  on 
the  order  of  ,5  seconds  with  approximately  1%  precision. 

The  authors  appreciate  the  assistance  rendered  by  the 
following  professors  in  various  phases  of  the  programs 
E.  A,  Milne  and  W.  W.  Hawes,  for  supervising  reactor  operations 
whenever  called;  J.  R.  Borstings,  for  advice  on  statistical 
matters;  and  H.  E.  Handler,  for  general  counseling  and 
assistance  throughout  the  program.  The  authors  also  wish  to 
acknowledge  the  assistance  Prof,  G,  W.  Rodeback,  thesis  advisor  $ 
gave  during  this  project. 
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I.  Introduction 

Comprehensive  investigations  of  isomeric  states  of 
isotopes  of  Hafnium  and  Ytterbium  have  been  reported  several 
times  since  the  discovery  of  nuclear  isomers  /l~5/o  The  results 
of  the  half-life  measurements  reported  by  E0  C0  Campbell  and 
P.  F.  Fettweis  in  1959  A/  are  summarized  in  the  following 
table,  along  with  the  corresponding  reactions  used  and  the 
gamma  energies  detected. 

reaction  half-life 


Ybl76  (n,tf)  Ybl77m  s  6.5  sec 
Hfl77  (n,tf)  Hfl73ra  =  3o5  sec 
Ta173  k  capture H£l7^m  -  4,8  sec 

Hf17^  (n,£)  Hf17^  -  18  06  sec 


gamma  energy 

.10  Mev 
.33,  .44  Mev 
.8  Mev 
.22  Mev 


Of  interest  is  the  apparent  difference  in  their  results  for 
the  half-life  of  Hf^7&ni^  although  this  may  have  been  within 
the  limits  of  the  experimental  errors,  which  were  not  reported. 
More  recent  work  (1962)  by  K,  F.  Alexander  and  H.  F0  Brinckmann 
/5,  indicates  half-lives  of  (4<>3  “  0el)  sec  and  (1896  t,  0o8)  sec 
for  the  Hfl77,178  (ir^Jf)  Hfl7&n,179m  reactions. 

The  objective  of  this  work  has  been  to  improve  the  precision 
of  the  half-life  measurements  of  Yb^-77m>  pj.£>178m^  and  Hfl79m0 
Measurements  were  made  by  scintillation  crystal  detection  of 
direct  gamma  decay  of  thermal  neutron  activated  samples.  The 
time  delay  between  activation  and  measurement  was  minimized 
through  the  use  of  a  rapid  pneumatic  transfer  system  which 
transfers  the  sample  from  the  activation  region  to  the  cotinting 
area.  Use  of  a  single  channel  analyser  and  fast  scaler  kept 
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dead  times  below  those  commonly  encountered  with  multichannel 
instruments 0  A  fast  continuous-- chart  recorder  was  employed 
for  accurate  recording  of  data0 

Data  was  analysed  in  a  computer  program  (FRANTIC  II,  by 
P„  Co  Rogers  of  M.I.To)  / 6/  which  performs  an  iterative  least 
squares  analysis  of  multiple  exponential  decays »  Several  runs 
have  been  made  for  each  nuclide,  not  only  to  check  reproducibility 
of  the  measurements,  but  also  to  determine  mean  half-life 
values  based  on  statistical  weighting  and  the  corresponding 
error  estimates  for  these  mean  half-lives 0 
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II.  Experimental  Equipment  (Schematic,  Fig.  1) 

A0  Transfer  System 

1.  The  rapid  pneumatic  transfer  system  was  designed 
and  built  for  use  in  the  AGN-201#  reactor  at  USNPGS  by 
Paas  and  Sullivan  /7/„  It  transports  a  sample,  which  is 
placed  in  a  small  light  container,  from  the  counting  area, 
about  eight  feet  from  the  reactor  shielding,  into  a  high 
neutron  flux  region  in  the  reactor  core.  Here  the  sample 
is  activated  by  slow  neutrons  and  then  propelled  under 
positive  pressure  back  into  the  counting  position 0  Transit 
time  from  core  to  counting  position  is  approximately  one~ 
tenth  second.  The  transfer  system  is  completely  closed, 
allowing  recycling  of  the  sample  without  disturbing  the 
detection  arrangement. 

2.  The  sample  container,  or  "rabbit ",  is  made  of 
balsa  wood  in  the  shape  of  a  bullet.  It  is  4° 5  cm  long, 
has  a  1.9  cm  diameter,  and  weighs  about  1.1  g  including 
the  sample.  The  isotopically  enriched  sample,  50  mg  in 
powdered  form,  was  placed  in  a  gelatin  medicinal  capsule 
which  was  then  secured  within  the  "rabbit".  Irradiation 
and  counting  of  a  dummy  "rabbit"  (less  sample)  indicated 
no  increase  over  background  rates. 

^The  AGN-201  is  a  graphite  and  polyethylene - 
moderated  research  reactor  capable  ''f  product ne  »  neutron 
flux  up  to  4  x  1010  n/cm^sec  at  the  AEC  license  limit  of 
1000  watts. 
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Figure  X • . 

Counting  System  and  Pneumatic  Transfer  System 


3.  Samples  were  procured  from  Oak  Ridge  National 
Laboratory  and  are  highly  enriched  in  the  nuclides  under 
investigation.  Isotopic  analyses  are  included  as 
Appendix  III, 

B.  Detector  System  (Fig.  2) 

1.  A  2"  x  2"  Nal(Tl)  crystal  was  used  as  the 
scintillation  detector.  The  crystal  used  was  a  Harshaw 
Model  SSS/2-I  integral  assembly  incorporating  a  Dumont 
6292  photomultiplier  tube  with  an  S-ll  response.  The 
assembly  is  rated  by  the  manufacturer  to  be  capable  of 
7.6%  resolution.*  Our  measurements  showed  the  resolution 
to  be  B.0%  at  the  recommended  1000  volts  anode-cathode 
operating  voltage.  During  counting,  the  crystal  was  seated 
flush  against  the  rabbit  tube  in  the  shielded  counting 
area.  The  shielding  material  was  lead  with  a  minimum 
thickness  of  two  inches.  The  counting  area  was  off-set 
about  16  inches  from  the  rabbit  tube  center-line  to  reduce 
the  background  rate  at  high  reactor  powers.  A  typical 
value  of  background  is  one  count  per  second  at  an  operating 
power  of  20  watts. 

2.  The  crystal-PM  tube  assembly  was  connected  to 
a  6AK5  cathode-follower  preamplifier  which  relayed  the 

*The  resolution  is  defined  as  the  ratio  of  the  energy 
width  at  half  maximum  to  the  mean  energy  of  the  .662  Mev 
Csl37  photopeak. 
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signal  through  50  feet  of  shielded  cable  to  the  input 
of  the  linear  amplifier* 

3.  High  voltage  for  the  photomultiplier  tube  was 
furnished  by  a  regulated  power  supply  with  a  specified 
stability  of  less  than  .0026$  output  change  per  line  volt 
input  change,  and  a  maximum  output  change  of  .2$  per  day. 

No  significant  H.  V.  drift  was  observed  during  the  decay 
runs.  Preamplifier  power  (*300v,)  was  supplied  by  a 
stabilized  power  supply  unit  in  the  single  channel 
analyser. 

C.  Single  Channel  Analyser 

A  commercial  combination  nonoverload  linear  amplifier 
and  integral-differential  discriminator  was  used* 
Continuous  amplification  (voltage  gain)  from  unity  to  960 
was  available  with  the  linear  amplifier.  Amplifier 
linearity  (section  III-A)  and  an  acceptable  signal-to-noise 
ratio  (45:1  for  Cs-*-37  photopeak)  were  verified  for  gain 
settings  in  the  range  used  for  all  runs.  The  discriminator 
unit  allowed  continuous  selection  of  pulse  heights  from 
0-100  volts  (normal  calibration  units)  and  differential 
pulse  selection  from  OtIO  volts.  The  whole  assembly  had  a 
resolving  time  of  approximately  one  microsecond. 

D.  Count  Recording  System  (Fig,  3) 

1.  The  output  of  the  discriminator,  a  20  volt  pulse, 
was  transmitted  to  a  decade  scaler  having  better  than  one 
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microsecond  resolving  time.  The  scaler  has  been  modified 
by  the  addition  of  a  6CB7  cathode  follower  which  acts  as 
an  impedance  matching  device .  Thus  when  the  scaler  was 
operating  in  a  continuous  mode,  scaled  output  pulses  could 
be  fed  to  the  chart  recorder  without  stopping  or  recycling 
the  scaler.  Scaled  output  pulses  could  be  selected  for 
each  2,  10,  102,  or  103  input  pulses  to  the  scaler 0  The 
scaled  pulses  were  sent  through  a  differentiating-clipping 
circuit  and  then  to  the  continuous  chart  recorder. 

2.  The  continuous  recorder  is  a  two»channel  chart 
type  with  variable  chart  speeds  of  1,  5*  25,  and  125  ram/sec. 
It  is  capable  of  recording  up  to  about  110  cps  with  a 
minimum  input  sensitivity  of  10  m.v.  For  this  application, 
output  pulses  from  the  scaler,  representing  a  scaled  number 
of  counts,  were  recorded  in  one  channel  and  a  calibrating 
time  mark  was  recorded  in  the  other. 

E.  Time  Calibration  System 

The  basic  instrument  for  this  system  is  a  crystal 
controlled  oscillator  designed  primarily  as  a  frequency 
standard  for  radio  frequency  calibration.  Frequency 
standards  of  10,  20,  or  100  kilocycles,  with  rated  stability 
of  *005%,  could  be  selected.  Pulses  of  desired  frequency 
were  transmitted  via  shielded  cable  to  a  five=decade ,  200 
kilocycle  scaler.  The  scaler  was  modified  to  act  as  a 
frequency  divider  by  tapping  off  the  third  or  fourth  stages, 
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thus  giving  a  wide  selection  of  time  calibrating  pulses 
(10"1,  1,  10,  etc,  seconds  per  pulse  interval),,  The  time 
pulses,  after  differentiation  and  clipping,  gave  a  single 
sharp  spike  per  scaled  number  of  input  cycles  and  an 
accurate  record  of  time  on  the  recorder  chart „ 


10 


III.  Calibration  and  Procedure 
A.  Energy  Calibration 

1.  A  planchette  containing  an  evaporated  Cs^V  sample 
was  used  as  a  calibrating  source.  The  known  $  energy 
(662  Kev)  and  long  half  life  (30  yre )  suited  it  to  this 
use.  The  planchette  was  secured  in  a  phenolic  holder  which 
was  machined  to  fit  over  the  crystal  case  (Fig.  4)°  The 
holder  and  crystal  case  were  scribed  so  that  a  reproducible 
position  was  available. 

2.  The  amplifier  linearity  was  adjusted  according 
to  the  manufacturer’s  instructions.  An  output  pulse  from 
a  mercury  switch  pulser  that  matched  the  pulse  shape  of 
the  signal  from  the  preamplifier  was  used.  At  the  same 
time,  for  convenience,  the  amplifier  was  calibrated  so 
that  the  lower  discriminator  (E)  settings  gave  the 
approximate  pulse  height  in  volts.  The  use  of  Cs^7  for 
an  energy  calibration  source  and  knowledge  that  the 
primary  gamma  energies  of  the  nuclides  of  interest  lay 
under  662  Kev,  led  to  the  choice  of  approximately  90  v. 

(.9  of  full  scale)  for  the  position  of  the  Cs^^  photoline. 
Gain  settings  on  the  amplifier  were  adjusted  so  that  the 
photoline  fell  at  roughly  90  v.  when  observed  on  a  cathode 
ray  oscilliscope.  Exact  position  of  the  photopeak  was 
determined  from  a  differential  energy  spectrum  and  found 
to  be  at  9S.5  volts  as  read  on  the  lower  level  discriminator 
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Calibration  Source  and  Holder 


scale.  Thus  the  discriminator  dial  setting  was  calibrated 
versus  gamma  energy.  A  check  on  the  linearity  of  the 
entire  system  and  verification  of  the  calibration 
described  above  was  obtained  by  measuring  the  photopeak 
position  of  the  364  kev  gamma  from  1^31 0  Results  are 


Fig.  5.  Energy  Calibration  Curve 

Periodic  checks  using  the  mercury  switch  pulser  as  the 
input,  and  an  oscilliscope  to  determine  pulse  heights  of 
the  amplifier  output,  indicated  no  departure  from  amplifier 
linearity  within  the  uncertainties  of  the  measurements 0 
Calibration  of  the  differential  discriminator,  or 
"window1' ,  gave  a  maximum  window  width  of  6.3  volts, 
corresponding  to  approximately  a  60  kev  energy  interval 
on  the  curve  shown  in  Fig.  5» 

3.  Daily  gain  adjustments  to  provide  energy  calibration 
were  made  by  a  "quick"  method  suggested  by  H0  E0  Handler  /6/„ 
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Initially,  this  method  required  the  selection  of  a 
reference  discrimination  setting  and  measurement  of  the 
corresponding  integral  rate  immediately  after  obtaining 
the  Csl37  spectrum,  A  value  near  the  midpoint  of  the 


(35,500  £  190)cpm 


dE 


yov  '  189, 5 V 

Discriminator  (E)  Setting 
Fig.  6.  Csl3?  Spectrum  Showing  Reference 
Voltage  and  Reference  Rate 


"valley"  between  the  Compton  edge  and  the  photopeak  was 
chosen  as  the  reference  discriminator  setting.  This 
choice  of  setting  makes  the  corresponding  integral  count 
least  susceptible  to  small  gain  changes.  At  this  setting, 
76  volts  on  the  lower  discriminator,  the  integral  rate 
(corresponding  to  the  area  under  the  curve)  was  measured 
as  (36,500  t  190)  counts  per  minute.  Once  this  rate  was 
determined  it  became  a  reference  for  future  measurements. 
Subsequent  energy  calibrations  were  made  by  first  adjusting 
the  amplifier  gain  to  give  the  approximate  desired  pulse 
height  of  the  reference  photopeak  as  measured  on  an 
oscilliscope.  Then  fine  gain  adjustments  were  made  at  the 
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reference  discriminator  setting  to  attain  the  reference 
counting  rate. 

B.  Gain  Stability 

The  fact  that  certain  photomultiplier  (PM)  tubes 
exhibit  marked  gain  shifts  (effectively  a  change  in 
multiplication)  at  high  counting  rates,  particularly  rates 
greater  than  5000  counts  per  sec  ,  has  been  reported  by 
several  investigators  /9,  10,  11/.  D.  F.  Coveil  and 
B.  A.  Euler  at  NRDL,  in  an  extensive  study  of  various 
brands  and  types  of  tubes,  found  that  of  the  tubes  shoxtfing 
gain  shifts,  all  required  an  hour  to  several  hours  to 
recover*  R.  Chery  corroborated  this  information  when  he 
found  PM  tubes  would  be  very  stable  for  pericds  of  about 
24  hours  if  they  were  first  exposed  for  20=30  hours  to 
high  (just  below  tube  dissipation  limits)  gamma  radiation 
levels. 

This  slow  tube  recovery  has  been  utilized  to  advantage 
in  this  experiment  by  presaturating  the  PM  tube  prior  to 
the  decay  runs.  The  sample  was  irradiated  at  the  same 
power  level  and  for  the  same  irradiation  time  as  for  the 
actual  decay  run.  Reference  rates  were  checked  after  this 
exposure  and  readjusted,  as  necessary.  After  decay  runs, 
reference  rate  checks  showed  less  than  1.5%  gain  shift 
(in  our  case  less  than  .75%  reference  rate  shift),  and 

*The  initial  gain  shift  occurred  almost  instantaneously. 
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tube  recovery  times  lasted  up  to  several  hours .  In  all 
cases,  recovery  times  were  very  much  greater  than  the 
maximum  counting  period  used  (five  minutes )c  The  effect¬ 
iveness  of  the  technique  was  substantiated  in  subsequent 
analyses  of  the  data  (Section  IV),  which  indicated  no 
appreciable  systematic  errors  as  would  have  been  the  case 
with  large  gain  shifts. 

An  alternate  method  for  checking  gain  shift  consisted 
of  placing  a  Cs^-37  sample  near  the  crystal  during 
measurement  runs  to  give  a  steady  high  background  rate. 
This  rate  (— -150  cps)  was  measured  before  and  after  each 
run,  at  operating  powers.  This  method  confirmed  that  gain 
shift  was  generally  less  than  1.0%,  and  data  were  not  used 
if  the  gain  shift  exceeded  1.5^»  Two  advantages  were 
gained  with  this  check.  The  first  was  that  the  measure¬ 
ments  could  be  taken  without  disturbing  the  crystal 
assembly.  Second,  this  allowed  a  visual  check  during  the 
run  by  superimposing  the  Cs-^7  pulse  over  the  measured 
pulses  on  the  oscilloscope.  In  this  manner  any  large 
short-term  shifts  during  measurement  could  be  observed, 
although  none  were  observed. 

Other  electronic  equipment  in  the  counting  system, 
once  stabilized  by  a  three  to  four  hour  warm-up  period, 
had  negligible  drift  over  the  periods  of  time  required  to 
make  one  or  several  runs. 
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Co  Measurement  of  Resolving  Time 

A  double  pulse  generator  was  used  to  determine  the 
resolving  time  of  the  counting  system,,  Initially,  each 
component  (amplifier,  discriminator ,  and  scaler)  was 
checked  separately  to  assure  specifications  were  meto 

A  test  signal  circuit  is  provided  in  the  preamplifier 
system  that  allowed  checking  the  system  resolving  time 
with  the  double  pulse  signal 0  The  results  of  measurements 
on  the  entire  system  gave  a  value  of  (1„4  micro-seconds 

for  the  resolving  time® 

Do  Activation  Parameters 

Activation  times  and  reactor  powers  were  varied 
according  to  the  half  life  measured  and  the  gamma  energy 
investigated o  An  effort  was  made  in  each  case  to  have 
an  initial  recorded  count  rate  of  5 » 000-10 9 000  counts  per 
second o  This  limitation  was  imposed  by  the  use  of  the  100 
counts  per  mark  scaling  factor  at  the  beginning  of  each 
run®  Rates  greater  than  10,000  counts/sec  approached  the 
recorder  frequency  limit  of  about  110  cycles/sec® 

The  Hafnium  activities  were  measured  at  gamma  energies 
of  (330  £  30*)  kev  for  Hf17gm  and  (215  £  30)  kev  for 
Hfl79mo  Hf17*m  axso  has  a  very  strong  (41$)  line  at 
213 o 4  kev  /ll /,  but  the  low  activation  cross-section  ratio 

jjc 

*  One-half  the  energy  "window",  rather  than  an  error 
estimate® 


1? 


/5 /  of  C'act(17Si/fr''  act (179)  s  °°54  for  our  sample 


precluded  the  use  of  this  energy  for  the  measurement 
of  the  shorter-lived  nuclide*  To  eliminate  possible 
degradation  of  the  results  for  the  longer  half-life  due 
to  pulse  pile-up*,  the  shorter-lived  nuclide  was  allowed 
to  decay  before  measuring  the  longer* 

In  general,  activation  times  of  approximately  one 
half-life  were  selected  in  order  to  obtain  at  least  50 fo 
of  the  saturation  activity  of  the  nuclide  of  interest, 
yet  reduce  contributions  to  the  counting  rate  from  any 
long-lived  impurities  present*  This  method  proved  inadequate 
for  the  investigation  of  the  Hf^Sm  nuclide*  Consequently, 
a  higher  reactor  power  for  a  shorter  irradiation  period 
was  chosen  to  give  a  ratio  of  initial  relative  activities 
of  at  least  3^1  for  the  Hff7$m  s  Hfl?9m  nuclides* 

The  integral  mode  of  operation  was  used  in  the 
Ybl77m  measurement*  The  low  cross-section,  about  7  barns 
for  Yb^^,  would  have  required  an  excessive  amount  of 
reactor  power  in  order  to  obtain  desired  counting  rates  in 
the  differential  mode*  The  integral  discriminator  (E  dial) 
was  set  to  eliminate  all  gamma  energies  below  about  95  kev* 

In  this  manner,  both  the  10/*.  and  212  kev  gammas  were 
measured  at  a  considerable  savings  in  reactor  power* 

*In  the  case  in  which  two  half-lives,  say  t^  and  9 
are  represented  in  the  source,  the  effective  pulse  intro¬ 
duces  three  false  half-lives  into  the  measured  activity: 


The  following  table  summarizes  the  activation 
parameters  and  the  gamma  energies  for  each  type  run 


Nuclide 

Hfl7&n 

Hfl79m 

Ybl77m 


Activation  Time 


Reactor  Power 


2  sec 
20 
5 


20  watts 

1 

10 


(330  Is  30)  kev 
(215  £  30) 
all  above  90 
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IV0  Data 


A0  Data  Reduction 

Data  tapes  consisted  of  two  channels  of  simultaneously 
recorded  information;  one  of  time*  the  other  of  scaled 
counting  marks  (Fig0  7)«  The  "counts"  channel  was  broken 
down  into  a  series  of  consecutive  time  intervals ,  starting 
with  intervals  that  were  short  compared  to  the  measured 
half -life  and  increasing  in  length  (time)  as  the  counting 
rate  decreased  and  approached  the  background  level »  Short 
intervals  were  chosen  so  that  the  sample  would  not  decay 
appreciably  over  any  given  counting  intervale  The 
approximate  expression  for  the  measured  rate* 

/ 1  f  ^ 


n  dt  * 


is  valid  only  for  the  case  of  At  small  compared  to  tj.; 
where  n  is  the  instantaneous  counting  rate  at  the  midpoint 
(t)  of  the  n  ^  interval  (At)  and  Rn  is  the  measured  number 
of  counts  per  interval  At.  In  each  of  our  decay  analyses* 
intervals  of  less  than  l/7  of  the  (shortest)  half-life 
were  used  initially 0  To  ensure  adequate  statistical 
accuracy  for  all  intervals  containing  more  than  2500  counts* 
25  marks  at  a  scaling  factor  of  100*  the  number  of  counts 
was'  estimated  to  the  nearest  50  counts 0  Similarly*  the 
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H 


FIGURE  7 

Data  Tape  (section  'of  tape  from  run  shovn  in  sample  output,  Appendix  l) 


number,  for  intervals  containing  625~2500  counts,  was 
estimated  to  the  nearest  25  counts «,  In  this  manner  the 
counts  in  each  interval  were  estimated  to  within  the 
standard  deviation  calculated  according  to  Poisson 
statisticso 

The  timing  error  was  estimated  by  measuring  the  drift 
of  the  calibrating  marks  in  the  "time”  channel  with  respect 
to  the  printed  grido  Typical  of  the  values  obtained  in 
these  measurements  (taken  over  the  entire  counting  period) 
was  1  x  10~3  se c/sec o 
Bo  Analysis  (least  squares  program) 

1 o  Theory 

The  data  for  the  two  component  decay  of  the  Hf 
sample  and  for  the  single  decay  of  the  Yb  sample  were 
analysed  with  the  aid  of  a  CDC  1604  computer 0  In  the 
iterative  least  squares  analysis  program  (FRANTIC)* 
the  input  data  points  (R^)  are  corrected  for  dead 
time  (T^)  and  background  (B)„  Then  each  point  is 
normalized  (XNORM)  as  desired  (in  this  case  XNORM 
s  loOO  was  used)  to  yield  the  correct  data  (An)  for 
the  data  point, 

(1)  A_  ,  Ob-  -  il  XNORM 

L1-^  -l 

A  statistical  weight  factor  (Wn)  is  then  applied, 

*This  discussion  will  necessarily  be  limited  to  the 
main  points  of  the  analysis „  For  further  clarification 
refer  to  item  #6,  Bibliography 0 
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which  is  of  the  forms 


1 


q£( total )n  (XN0RM)2 


where  &?( total )ns(  the  total  variance  for  the  pointy 


consists  of  the  sum  of  the  variances  for  1)  count  rate 
uncertainty  2)  background  uncertainty  3)  uncertainty  in 
dead  time  and  4)  uncertainty  in  the  time  length  of  the 
counting  interval » 

For  finite  observation  intervals  (DTn)s,  the  expression 
for  the  activity  (ACn)  resulting  from  N  radioactive  decays 
is : 


This  exact  expression  is  not  properly  evaluated  by 
the  computer  for  small  values  of  >|DTn  so  that  equation  (3) 
is  expressed  in  terms  of  a  hyperbolic  sine  function  and 
expanded  to  give  the  expression: 


(4)  ACn  s  £  Ao„  e 

»  g  i  i 


In  the  least  squares  fit  analysis  this  expression  is 
used  for  values  of  X*  DTn  up  to  lo0  and  the  original 
expression  (3)  is  used  for  values  of  ^;DTn  greater  than 
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l„0o  The  X-,  values  are  then  held  fixed  and  the  above 
expressions  are  linear  in  the  An  coefficients 0 

In  the  iterative  least  squares  analysis  the  exact 
expression  (3)  is  expanded  in  a  Taylor  series  about  the 
point  defined  by  the  values  of  A0«^and  X-  f rora  the  preceding 
iteration 0 

The  iterative  process  is  continued  until  all  parameters ? 
as  well  as  the  weighted  variance  of  fit  (the  square  of  the 
standard  deviation  of  the  distribution  of  residuals  about 
zero),  deviate  from  their  value  in  the  preceding  iteration 
by  not  more  than  one  part  in  10^,  or  until  25  iterations 
have  been  completed 0 

A  sample  of  the  computer  output  from  a  measurement  of 
Hfl?9m  t[S  included  as  Appendix  IQ 
2«  Goodness  of  Fit 

Of  particular  interest  in  this  program  is  the  goodness 
of  fit  information  included  in  the  printout.  the 

weighted  variance  of  fit  (VAR),  and  a  histogram  of  the 
distribution  of  the  residuals  (of  the  experimental  curve 
about  the  calculated  curve)  are  given.  In  using  the  usual 
^  2  VSo  degrees  of  freedom  tables,  the  weighted  variance 
of  fit  (VAR)  times  the  degrees  of  freedom  (DF)V  will  give 
a  more  realistic  measure  than  X  ^  of  the  probability  of 

*In  this  case  degrees  of  freedom  is  defined  as  the 
number  of  input  data  points  minus  the  number  of  parameters 
allowed  to  float. 
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performing  a  better  experiment  or  finding  a  better  fit 
with  different  parameters  or  a  different  hypothesise  This 
is  true  since  the  uncertainties  in  dead  time,  background, 
and  length  of  time  interval  are  considered  in  the 
calculation  of  VAR,  whereas  they  are  not  in  calculating  X. 

As  an  example,  using  the  output  data  shown  in  Appendix 
I,  VAR  x  DF  =  27.55#  (with  29  degrees  of  freedom).  This 
gives  a  91.5$  "assurance”  that  the  calculated  fit  is  unique 
for  the  given  data,  i.e0,  that  the  hypothesis  proposed  is 
correct.  The  histogram  of  the  distribution  of  the  residuals 
is  approximately  Gaussian,  as  it  should  be.  In  this  case, 
had  a  twro  component  decay  mode  been  proposed  in  the  initial 
estimates,  the  anlysis  would  have  either  rejected  the 
second  mode  or  given  a  large  VAR  indicating  a  false 
hypothesis. 

In  general,  an  arbitrary  criterion  for  rejecting  data 
was  established.  A  significance  level  (OC )  of  .10  was 
chosen  for  the  individual  runs.  This  means  that  we 
accepted  a  10$  chance  of  making  an  error  when  rejecting  a 
result  as  false;  or,  conversely,  only  10$  of  the  time, 
over  a  series  of  many  runs,  would  true  results  fall  outside 
of  our  criterion.  In  the  example  mentioned  above  the 
value  of  (X  for  this  run  would  be  .915,  meaning  a  91.5$ 
chance  of  being  incorrect  if  the  results  were  rejected 
as  false.  As  an  example  of  a  run  that  was  rejected;  the 
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VAR  x  DF  was  61.65.  From  the  X  2  vs.  degrees  of  freedom 
tables,  we  find  that  there  is  only  a  5$  chance  of  being 
incorrect  in  rejecting  this  result  /l4/o 
3 .  Method  of  Weighting  Results 

Results  of  the  half  life  measurements  were  weighted 
according  to  the  "goodness  of  fit"  of  the  least  squares 
analysis . 

The  weight  factor,  ,  applied  to  each  result  is: 


1 

where:  s  (VAR  x  AM”1)^ 

VAR  -  weighted  variance  of  fit 
AM  ”1  s  the  inverse  of  the  least  squares 
matrix  in  the  analysis. 

This  value  is  the  value  given  as  the  half  life 
"sigma"  (App.I)  in  the  computer  program  output  for  the 
particular  run  considered. 

Using  this  weight  factof  the  half  life,  ti,  is  found 


by: 


where : 


n 


the  half  life  for  the  i^*1  run 
number  of  runs  used. 
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The  standard  deviation  of  the  half  life,  (A(ti),  is 


calculated  from: 


CMti) 


n 

.1 

Is  1 
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V.  Results 


In  the  analysis  of  the  data  the  mode  of  decay  (i0e0,  two 
component  or  one  component)  used  in  the  program  depended  on  the 
half  life  being  measured.  In  the  case  of  Hf^75m  a  two 
component  decay  was  assumed.  For  Hf^79m  and  Yb^77m  single 
component  decays  were  used.  This  required,  in  the  case  of 
Hfl79m}  that  the  shorter  half  life  be  allowed  to  decay  to  a 
small  value  (less  than  1%  of  total)  before  measurements  were 
taken. 


The  tabulated  results  of  the  individual  acceptable  runs 

for  the  measurements  of  the  half-lives  of  the  isotopes  under 

consideration  are  as  follows; 

Hfi77  (n,  8s)  Hf17^ra  (  .326  Mev  transition  measured) 

Half-life 

3.774  sec. 

3.996 
3.952 
3.917 


,177  sec: 

,209 

.057 

,060 


Hf17&  (n,8M  Hf179m  (# 217  Mev  transition  measured ) 
Half-life 

15,455  sec. 

15.519 
15.509 
15.574 
15.373 

Yb^-76  (n,  ^  )  Yb**-77m  (0212  Mev  transition  measured) 

Half-life 


6.564  sec. 
6.470 

6.466 

6.542 


,092  sec, 

>073 

,165 

>099 

,090 


,092  sec, 

,046 

,046 

,091 


2g 


The  weighted  averages  resulting  from  these  measurements  are; 


Isotope  Half-life 

(3o92  t  o05)seco 

(13o49  t  e04)seco 
(6o49  t  o03)seco 

The  value  of  IS ,49  sec  for  the  half-life  of  Hf^^^m  agrees 

with  the  value  of  1S»6  t  „£  reported  by  Alexander  and 

Brinckmann „ 

The  value  of  3° 92  sec»  for  the  half-life  of  Hf^Sm 
falls  almost  exactly  between  the  measurements  of  3°5  sec0 
reported  by  Campbell  and  Fettweis,  and  the  4<=3  sec0  value 
reported  by  Alexander  and  Brinckmann „ 

The  value  of  6 <,49  sec,  for  the  half-life  of  Yb^^^m 
is  in  good  agreement  with  the  value  of  6 0  5  sec,  also  reported 
by  Campbell  and  Fettweis „ 

Assuming  that  there  are  no  appreciable  systematic 
errors,  it  would  appear  that  the  present  system  of 
recording  and  analysing  decay  data  is  capable  of  good 
precision,,  Thus  we  are  led  to  prefer  the  single  channel 
device  for  the  following  reasons ; 

1)  The  system  dead  time  is  considerably  less  than  that 
encountered  with  multichannel  analysers  (1„4  micro-sec  vs„ 
20-50  micro-sec  for  many  multi-channel  devices )0  Smaller 
dead  time  allows  working  with  higher  activities,  hence 
better  precision  in  counting  data0 


Hfl7$m 

Hfl79m 

yb177m 
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2)  The  dead  time  factor  and  overall  timing  error  can 
be  accurately  determined  for  the  single  channel  instrument 0 
Variations  which  occur  in  "live"  (on)  times  and  switching 
(off)  times  between  channels  make  evaluation  of  actual 
counting  times  and  dead  times  a  difficult  task  with  multi¬ 
channel  devices. 

It  is  felt  that  the  present  system ,  with  the  FRANTIC 
computer  program,  can  be  effectively  used  to  measure  half- 
lives  on  the  order  of  0«5  seconds. 
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ANALYSIS  OF  HFU79M)  3  MAY  63  RUN  NO*  2  lW,  16  SEC 
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APFENDI5C  I 


Appendix  II 

Isotope  Energy  Levels  and  Decay  Schemes 
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APPENDIX  III 


Isotopic  Analysis  of  Samples 

I,  Analysis  of  Hafnium  Sample 
(as  reported  by  Oak  Ridge  Analysis) 


Element 

Hafnium 

Isotopic  Analysis 

Isotope 

177 

Isotope 

Atomic 

Precision 

percent 

Series 

G.  P. 

174 

trace  only 

Sample 

327  (a) 

176 

<>99 

♦ 

W 

o07 

177 

59,08 

T> 

o  52 

173 

31,73 

o  49 

197 

4,53 

«■ 

,17 

130 

3  c  68 

+ 

oil 

The 

limits  quoted 

above  are 

an  expression  of 

the 

precision  of  this  measurement  only.  The  error  is  estimated 


at  less  than  1 %  from  known  sources  of  systematic  errors. 

Spectroscopic  Analysis 


Element 

Percent 

Cu 

Ool 

Fe 

0,2 

2.  Analysis  of  Ytterbium  Sample 
(as  reported  by  Oak  Ridge  Analysis  ) 
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Element  Ytterbium 


Isotopic  Analysis 


Isotope 

176 

Isotope 

Atomic 

percent 

Precision 

Series 

IU 

Sample 

1061  (a) 

168 

A 

o 

o 

170 

A 

o 

O 

171 

oil 

172 

026 

i  o05 

173 

o 

o 

1  o05 

174 

1.8l 

t  o05 

176 

97o57 

i  .i 

The  limits  quoted  above  are  an  expression  of  the  precision 
of  this  measurement  only*  The  error  is  estimated  at  less  than 
1%  from  known  sources  of  systematic  errors . 

Spectroscopic  Analysis 


A1 

o 

o 

V 

Ni 

4r\ 

O 

o 

V 

Gd 

< 

o  004 

Ba 

<  .02 

Pb 

<  .1 

Tb 

< 

.05 

Be 

<  .001 

Si 

<  .05 

Dy 

< 

ol 

Ca 

<  .07 

Sn 

C  .05 

Ho 

< 

o02 

Co 

<  .05 

Ti 

<  .02 

Fr 

< 

o004 

Cr 

A 

o 

o 

V 

<  .02 

Tm 

< 

.02 

Cu 

UA 

O 

• 

V 

Zr 

<  .1 

Lu 

< 

.018 

Fe 

<  .02 

Y 

<  o  004 

K 

<  .01 

<~  No  spectrum  line 

La 

<  o01 

Li 

<  .01 

visible 

o  Probably 

Ce 

<  .2 

<  ol 

absent, 

definitely 

Mg 

<  .02 

Pr 

Mn 

<  .02 

less  than 

value  given 

Nd 

<  ol 

Mo 

c  .05 

Sm 

<  a 

Na 

<»05 

Eu 

<  .004 
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